Spectral functions in the sigma-channel near the critical end point 
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Spectral functions in the a-channel are investigated near the chiral critical end point (CEP), that 
is, the point where the chiral phase transition ceases to be first-ordered in the {fj,, T)-plane of the 
QCD phcise diagram. At that point the a meson becomes massless in spite of explicit breaking of 
the chiral symmetry. It is expected that experimental signatures peculiar to CEP can be observed 
through spectral changes in the presence of abnormally light a mesons. As a candidate, the invariant- 
mass spectrum for diphoton emission is estimated with the chiral quark model incorporated. The 
results show the characteristic shape with a peak in the low energy region, which may serve as a 
signal for CEP. However, we find that the diphoton multiplicity is highly suppressed by infrared 
behaviors of the a meson. Experimentally, in such a low energy region below the threshold of two 
"^i^" ' pions, photons from tt" 27 are major sources of the background for the signal. 
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I. INTRODUCTION 



The spontaneous breaking of the chiral symmetry has fundamental significance in understanding the non- 
perturbative nature of hadron dynamics. It has been argued that the chiral symmetry can be restored at sufficiently 
' high temperature and/or high baryon density by means of effective models and lattice calculations (see Ref. |0| for a 
f — , state-of-the-art review of lattice results). In QCD with two massless (massive) quark- flavors at zero baryon density, 
psj ■ the chiral phase transition is supposed to be a second-ordered one (crossover) according to the universality argu- 
0^ ', ment [||. Once an adequate chemical potential /i for the baryon density is introduced, the chiral phase transition can 
■ become first-ordered for small quark masses, as suggested by several model studies ^, ||, ^, 01- 
' If the strange-quark mass is large enough to make the chiral restoration at zero baryon density a continuous 
transition, as found in the lattice calculation with staggered fermions the first-ordered line should terminate at 
' some point in the (/i, r)-plane of the QCD phase diagram. This terminal point is called the chiral critical end point 
Oh, (abbreviated to CEP in this paper). Two minima of the effective potential degenerate right at this point, so that 
O ^ the curvature around the minima vanishes. This results in the appearance of the a meson with zero screening mass, 
D , even though the pions are still massive due to explicit breaking of the chiral symmetry. This is the reason why much 
• attention has been paid to physical consequences around CEP from not only the theoretical but also the experimental 
^ [ point of view @, ||, |[ |l^, |ll[ As a matter of fact, future experiments of the heavy-ion-collision planned in GSI 
• ^ • with 30 A GeV energy and in JHF @ with - 25 A GeV energy will explore a lot about the high density nature of 
QCD, including physics around CEP. 

It would be interesting to consider the possibility to detect such light a mesons directly in heavy-ion-collision 
experiments. The a meson in a hot medium goes through such processes as cr ^ 27, a 27r, a 2(T, a — > NN, etc. 
It is expected that the measurement of two photons (diphotons) can tell us prosperous information on the transient 
thermal medium because electro-magnetic probes, such as leptons and photons, hardly receives rescattering in the 
relatively small system produced in a collision. Formerly the diphoton measurement was proposed as a candidate of 
the implements to see the spectral changes near the chiral restoration: The a mass is so reduced around the chiral 
transition temperature that the a meson cannot decay into two pions and thus the spectral function in the cj-channel 
must be significantly enhanced |l^, |l^, |2^, |2^, |2^ . In the present paper we apply the idea to the spectral 

changes near CEP. Actually strong vestiges of CEP can be anticipated since the spectral function has not just an 
enhanced peak but a pole contribution then. We are going to estimate the diphoton spectrum taking advantage of the 
universality argument around CEP. However, the information derived from the universality is not sufficient to obtain 
the diphoton emission rate. For that purpose, we must take account of each dynamical processes in order to acquire 
the spectral function. 

In the construction of the spectral function we would emphasize the following viewpoint: The cr mass and the scalar 
meson condensate are not sensitive to the detail of dynamical processes. Their smooth decreases at finite temperature 
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FIG. 1; Schematic picture of the QCD phase diagram. CEP is an abbreviation for the critical end point, which is located at 
(/iE,7E). The dotted line with an arrow indicates the evolving direction considered in the present analysis. 



can be provided by almost any thermal fluctuation. Therefore we can parametrize the behaviors a priori resorting to 
the universality argument or knowledges inferred from model studies, as will be performed in Sec. The width and 
the amplitude (residue of the propagator), on the other hand, strongly depend on dynamical processes, or in other 
words, depend on what kind of loops are taken into account. Thus we must evaluate them by considering respective 



processes, as will be calculated in Sec. [11 



Once we have the spectral functions in the cr-channel, we can investigate the effects on the diphoton emission 
coming from the spectral changes near CEP. The dominant contribution from the a pole brings about a characteristic 
shape in the diphoton spectrum, which reflects how long the system stays around CEP. In contradiction to a naive 
expectation as mentioned above, we will find in Sec. IV that the a pole contribution is strongly suppressed by the 
small residue in infrared regions. Thus the diphoton multiplicity is weaken due to infrared behaviors of the a meson. 
Finally we summarize our results and discuss the possibility to detect CEP via the diphoton measurement in Sec. 



II. PARAMETRIZATION NEAR THE CRITICAL END POINT 



In the vicinity of CEP located at (/ig, Tg), the critical behaviors are predominantly described by the lightest mode, 
i.e. the cr meson (long-ranged correlation in the scalar- isoscalar-channel) alone. Since the cr meson is a scalar particle 
associated with the Z(2) symmetry, the universality argument tells us that the chiral phase transition belongs to the 
same universality class as that of the 3-d Ising model. It is worth noting that the chiral restoration at CEP would 
become a genuine phase transition, namely, a second-ordered phase transition, even though it is a crossover at lower 
baryon chemical potential. Then singularities around that critical point are characterized by the critical exponents as 
functions of the distance from the critical point. In general this distance I in the (/^,r)-plane becomes an admixture 
of both temperature-like (denoted by t) and magnetic-like (denoted by h) variations in the Ising model: The Z(2) 
symmetry is only exact right at CEP and the effective potential is distorted away from that point. In the QCD phase 
diagram the Hike direction, along which deformations of the effective potential should preserve the Z(2) symmetry, 
is tangential to the first-ordered line, as depicted in Fig. [|. 

The critical exponents for the 3-d Ising model have been well examined analytically as well as numerically p3| . The 
inverse correlation function and the spontaneous magnetization would vanish as ~ ^0.63 g^^^ ^ ^ t^'^^ 

along the t-direction and as ~ h"^^^^ ~ ft.^-^" and m ~ h^^^ ~ along the /i-direction, respectively. It turns 
out from these numerical values that the magnetic singularities should dominate as long as I contains a non- vanishing 
component in the /i-direction, or in other words, they should dominate unless the direction specified by I is parallel 
to the pure t-direction tangential to the first-ordered line. As a result, singularities in almost all directions can be 
effectively described by the magnetic exponents, as discussed in Ref. Then the mass of the a meson, which can 
be interpreted as the inverse correlation length in the Ising model, vanishes near CEP as - /"'/'^^ while the scalar 
condensate drops off as ~ l^^^ as though it would be a spontaneous magnetization. The pion mass TOtt, on the 
other hand, can be approximately regarded as constant up to the critical point, which is common in several model 
studies 0, m. 

For simplicity we restrict our discussion throughout this paper only to the case with the chemical potential fixed 
at /X = /lE (the direction indicated by the dotted line in Fig. |l|). In this case, I is proportional to T — Te, and we can 
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parametrize the (screening) masses simply as follows: 



T 



m,(r)=™;, m,(T) = r<|l- j I , (1) 

with h'/ (3S — 0.403. Also the scalar condensate, that is, the counterpart of the spontaneous magnetization in the Ising 
model can be parametrized as 

■ \ 2^1 1/5 



MT)^vo + viT) = af: + {l-a)f:^l- i^—) j , (2) 

with 1/(5 = 0.208. Here vq is the condensate remaining at CEP due to explicit breaking of the chiral symmetry. It 
should be noted that the scalar condensate can be regarded as the pion decay constant under some approximations. 
Accordingly we denote the condensate as /ttCT) in Eq. (||). Indeed, this identification for the pion decay constant is 
used in the choice of the parameter sets, as stated below. 

The universality argument tells nothing about the values of m*, m*, /* (the pion mass, the a mass, and the 
pion decay constant at (fi — fiE, T — 0) respectively), and vq. In the present analysis these are all treated as free 
parameters, which should be arranged by hand or by using some models. Here, avoiding artifacts inherent in any 
model study, we simply employ the typical parameter sets as follows, that is, we take two extreme cases, one with no 
change of the ^-quantities and another with a relatively large change of the ^-quantities: 

(CASE-I) m; = 140 MeV, = 600 MeV, f* = 93 MeV, 

(CASE-II) ml = 140 MeV, = 300 MeV, f* = 46.5 MeV. (3) 

As for the choice of a, for the moment we will show the results of the spectral functions and the diphoton emission rates 
only for a specific choice of a = 1/2; the condensate at CEP decreases up to the half of its value at zero temperature. 
Finally we will present in Fig. ^ the results for a variety of a ranging from 0.2 to 0.9. The qualitative consequences 
are hardly amended by a change of a, though the absolute amount of the diphoton yield depends on a. 

We should remark that the choice of CASE-II may have a relation to the Brown-Rho scaling hypothesis in dense 
matter ||2^ , i.e. m*^/miy ~ /^//tt — 'i'(PE) where pE is the baryon number density corresponding to the chemical 
potential /^e- CASE-II implies the choice of $(yOE) = 1/2. Taking into account the fact that the partial restoration of 
the chiral symmetry might be experimentally observed in nuclei [p5| , p6| , we can anticipate that such reduction up to 
the half may happen in the vicinity of CEP. 



III. SPECTRAL FUNCTIONS NEAR THE CRITICAL END POINT 

From the experimental point of view, functional forms of the meson masses and the condensate assigned by Eqs. (|l|) 
and (^ are not sufficient to describe physical decay processes, that is, the decay width and the amplitude must be 
evaluated. These informations are contained in the spectral functions. Actually the spectral functions in the cr-channel 
are necessary to estimate the diphoton emission rates. Therefore, in this section, we are going to draw our knowledge 
on the spectral functions in the a-channel, limiting our consideration to the case with the spatial momentum fixed at 
zero. The spectral function in the cr-channel is defined as 

1 1 ImH^ (p) 

Pa{p) = ^-Iv^D^ip) - -- • (p2 _^2_ Rcn^V + (ImHR(p))2 ■ ^"^^ 

D^{p) and n^(p) are the retarded Green's function and the retarded self-energy in the cr-channel, respectively. In 
principle, all we have to do is calculate the self-energy n^(p) near the critical point as performed, for example, in 
Ref. 

In place of performing some resummation procedures, we calculate only the imaginary-part of the self-energy using 
the propagators with the parametrized masses. One-loop diagrams shown in Fig. ^ describe the decay processes 
a 2a and a — > 27r. We would emphasize that this is not a simple perturbative expansion but rather a resummed 
calculation since the propagator incorporates the interaction in the mean-field of masses given by m7r(r) and m^(T). 
Also the full-vertices represented by boxes in Fig. || can be inferred from the effective potential (Appendix ^ for detail 
arguments). 
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FIG. 2: Diagrams for the self-energy with non-vanishing imaginary-part up to the one-loop order. Boxes stand for the full- 
vertices inferred from the effective potential. 



Here we limit ourselves to the case in the back-to-back kinematics, that is, diphoton emission with zero spatial 
momentum. Then we obtain 



/ 4m2 



Imn?(cj, 0) = - 1 f 1 + 2nB{uj/2) e{uj - 2m,) 



327r V uj 



3t/2 4^2 



(TTTTT 



327r 



(l + 2nB(t^/2)) 0(t^-2m^), (5) 



where nB(w/2) stands for the usual Bose distribution function and 9{lj — 2m) is the step function. The expressions 
for the vertices Vaaa and V^jr^ are given in terms of m.„{T), m^lT), and f-n{T) in Eq. (A_3) in Appendix 



We have neglected processes involving quark (nuclcon) loops such as cr ^ qq (a ^ NN) because of physical and 
technical reasons. First of all, it should be negligible at CEP since the constituent quark is still heavy (~ 100 MeV), 
while the cr meson becomes light. Furthermore, the inclusion of quarks needs the wave- function renormalization, in 
which we cannot avoid technical complications and ambiguities. 

It is worth noting that the parametrized mass is the curvature of the effective potential, i.e. the screening mass. 
The pole mass appearing in the calculation of the self-energy might be different from the parametrized one due to 
the Lorentz anisotropy at finite temperature and/or density. In Ref. [ p7| it is clarified that the anisotropy results in 
the dynamical critical exponent z ~ 1 -|- 1/27 for the (/)^-theory slightly above Tc at zero density. The underlying 
physical contents of massless a mesons near CEP are effectively the same as the critical phenomena described by the 
^■^-theory. Thus we extrapolate the results of Ref. |^ and suppose that the effect of the dynamical critical exponent 
is small enough to be negligible even below Tc at finite density. This presumption of neglecting the dynamical critical 
exponent formally corresponds to a specific choice of renormalization conditions like the fastest apparent convergence 
(FAC) condition adopted in Hatree-type resummations [ p8[ , in which the coupling constant contained in models should 
depend on the temperature, in principle. 

Also wc briefly comment upon the results of Ref. |^ . The authors argued that the pole mass of the a meson within 
the NJL-model remains finite even right at CEP in the leading order of the 1/Nc expansion. This is not because of 
the treatment of the zero-point energy as discussed by the authors, but because of the lack of infrared dynamics. We 
have found that the massless a loops dependent on the external momentum are essential ingredients to describe the 
critical properties psj . 

Now that we have the imaginary-part of the self-energy, we can construct the real-part from its imaginary-part via 
the dispersion relation, i.e. 

Ren^(w, 0) = -V r As I"in^(Vf'0) ^ (subtractions), (6) 

apart from subtraction factors which are to be determined by the renormalization conditions, in principle. V stands 
for the prescription of Cauchy's principal value. As for the subtractions, meson one-loops need only one subtraction, 
which is absorbed in the mass-renormalization. We can fix the subtraction factor without ambiguity by demanding 
the a mass to be given by Eq. (jl|), neglecting the effect of the dynamical critical exponent. If quark loops were 
taken into account, they would need two subtractions in the one-loop order, which are absorbed not only in the mass 
renormalization but in the wave-function renormalization also. This makes the actual computation quite intricate, 
though the quantitative behaviors are essentially described by meson loops only. 

Then we can evaluate the spectral functions in the cr-channel by substituting ImnJ^(u;) and Ren5(ct;) for the 
definition of Eq. (^. m? in Eq. (^) is irrelevant since it can be absorbed in the subtraction factor. Throughout 
this analysis we set the location of CEP at (u^ = 700 MeV, Te = 140 MeV) as suggested by Refs. [|, |l|. The 
resulting spectral functions are shown in Fig. for the cases of T = 130 MeV and T — 139 MeV. It is clear from the 
shown spectral functions that the dominant contributions in the low energy region arise from the hght a pole and the 
continuum for the decay process a 2a. The threshold of the continuum contribution lies in the point uj = 2m^(T). 
We note that the process tr ^ 27r has only slight effects on the spectral function at T = 139 MeV simply because the 
full- vertex Vo-tttt becomes decreasing when the a mass approaches the pion mass, as can be seen in Eq. ([A3|). 
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FIG. 3: Spectral functions in the a-channel near CEP located at (^e = 700 MeV,TE = 140 MeV). The left figure is for the 
parameter set of CASE-I ($(pe) = 1) and the right for that of CASE-II ($(pe) = 1/2). 
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FIG. 4: Strength of the pole contributions. 



To investigate the spectral enhancement in a more quantitative sense, we define the strength of the pole contribution 



as 



(7) 



where is the temperature at which ma-{T = T^^) ~ u is satisfied. Then the pole part of the spectral function can 
be written as = Zc{uj)5{lj — m„{T)). In Fig. ^ we plot the strength of the pole contribution as a function 

of LJ. As is clear from the figure, the strength becomes vanishing at = 0. This property is understood immediately 
from the dispersion integral of Eq. (^. The differentiation of the integral diverges infraredly (s ~ 0) as a; ^ 0. As a 
result of the infraredly diverging denominator, the spectral function becomes zero at w = 0, which is consistent with 
the general property of the spectral function p{ljj) = — p(u;), for mesons.^ We will discuss this behavior of the pole 



strength in connection with the diphoton observation later in Sec. [V 



^ The author thanks T. Hatsuda for pointing it out. 
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FIG. 5: Diphoton emission rates for 6 — 0. The left figure is for the parameter set of CASE-I ($(pe) = 1) and the right for 
that of CASE-II ($(pb) = 1/2). 



IV. THE DIPHOTON EMISSION RATE 



One of the most promising candidates to see the spectral changes in the cr-channel is the invariant-mass spectrum 
of diphoton emission from the a meson. The emission rate per unit space-time volume is given by (Appendix ^ for 
the derivation) [p9| 



dN 



g/3cj cosh 2 ' 



(8) 



with the effective coupling ga^yj of the decay process a 27. Although the expression is proportional to uj'^, the 
singularities in g„^j(uj) just compensate for in the presence of a thermal medium and uj'^ga-yy{w) takes a finite 
value in the limit of w — *■ 0. The actual evaluation of (70-77 (w) as a function of lu is illustrated in Appendix |^. 6 is the 
fluid rapidity related to the fluid velocity v hy 9 — arctanhu arising from the Lorentz boost. 

Since they are proportional to each others, the gross features of the multiplicity are the same as those of the spectral 
function. The dominant contribution, in fact, comes from the a pole when the system lies in a state close to CEP. 
The results are presented in Fig. || for 6* = 0, in accord with the spectral functions in Fig. ||, respectively. 

We must take account of the space-time evolution of the system for the calculations to be compared with exper- 
imental outputs. Here we adopt the simplest situation for the space-time history, though it is not quantitatively 
realistic but qualitatively acceptable. 

First of all, we fix the baryon chemical potential at fi — fi^ during the evolution until freeze-out. This approximation 
makes the entropy per baryon deviate at most ^ 10% depending on the nucleon mass around CEP when the equation 
of state for the ideal gas is used. Considering that neither the nucleon mass nor the equation of state is precisely 
known around CEP, this approximate path with /i — ue is the best we can do for the qualitative study. 

Secondly, following the analysis presented in Ref. ||29|], we presume that the hydrodynamic evolution would obey 
the Bjorken's scaling solution In a quantitative sense the transverse expansion should be taken into account 

and moreover the boost invariance might be unsatisfied under the experimental conditions where CEP is concerned. 
Nevertheless, the analysis based on the scaling solution will provide us with the qualitative pattern of the diphoton 
spectrum when CEP is passed through. 

Then the invariant-mass spectrum with the transverse mass Aft, the transverse momentum pt and the momentum 
rapidity Y of diphotons fixed is given by 



1 



dN 



(Area) dM^d^ptdY 



Pt=Y=0 



dr T / dy 



dN 



d'^xd'^p 



(9) 



p=(M,()) 



where (Area) is the constant transverse area of the hadronic gas and y is the space-time rapidity which is equal to the 
fluid rapidity 9 in Bjorken's scaling solution, rj and Tf are the initial and the final time of the evolution; ti is taken as 
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FIG. 6: The temperature as a function of time derived from the ansatz given by Eq. (^3|) with the initial condition Ti = Tb = 
140 MeV at Ti = 1 fm. 



1 fm as usual and T{ is defined as the time when the system is cooled down below the freeze-out temperature. We fix 
the freeze-out temperature around CEP as T{ ~ 120 MeV according to Ref. In order to accomplish the integration 
with respect to r, we need the temperature at given t. That is determined by the following argument, as discussed 
in Ref. ||3^. Bjorken's scaling solution leads to the entropy density as a function of time as 

44 = ^, (10) 

S[Ti) T 

from which we can read the temperature once we have a relation connecting the entropy density to the system 
temperature. For the ideal gases, the entropy density is given by 

.H(r) ^ 3 . 14^^ + 4 • ^{t/ ^ (ln[l + e-«^-^-)] + ln[l + e-^(^^+--)] ) }, (11) 

for the hadronic gas composed of pions tt^, tt^, and nucleons p, n. Pions can be treated as massless particles, while 
Ek ^ ^/WTm^ is the energy of nucleons with, the mass in a medium. For the quark-gluon plasma composed of 
massless gluons and quarks with two flavors d at quark chemical potential ^e/3, the entropy density is given by 



In the intermediate temperature we can smoothly interpolate between them as follows: 

= li^QiT) + sh{T)) + i(sQ(r) - sh{T)) tanh(^^^^^ , (13) 

where AT controls the strength of the (deconfinement) transition. It is worth noting that the change of the entropy 
density is mainly attributed to the liberation of the color degrees of freedom and thus it has little to do with the 
chiral phase transition, at least in principle. There remain subtleties in this respect because no well-defined indicator 
of confinement is established so far If the deconfinement transition is really related to the chiral dynamics, as 
often said, the change of the entropy density at CEP {i.e. the terminal point of the first- ordered transition) could be 
almost discontinuous, that is, AT is small or zero in effect. Here we set AT — 1 MeV. 

The temperature as a function of time is shown in Fig. ^ with the initial condition T = Te = 140 MeV at ri = 1 fm. 
We take this initial condition, in which the gas is equilibrated right at CEP, for we have parametrized the masses 
below Te and have no idea about the masses above CEP. The dotted curve stands for the temperature evolution in 
the case when the nucleons are infinitely heavy. In CASE-I and CASE-II the nucleon mass M* is chosen in accord 
with the Brown-Rho scahng hypothesis, that is, M* ~ 1 GeV for CASE-I and M* ~ 0.5 GeV for CASE-II. 
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FIG. 7: The multiplicity of diphoton emission near CEP. The left figure is for the parameter set of CASE-I ($(pe) = 1) and 
the right for that of CASE-II ("I>(/5e) — 1/2). The dotted line indicates the position of the pion mass m„o = 135 MeV. 



We note that the ansatz given by Eq. (|_3|) is impUcitly based on a specific choice of the initial entropy density 
s(Te) = (sH(TE) + SQ(rE))/2. The gentle slope of temperature until r ~ 3fm (1.2fm) in CASE-I (CASE-II) is caused 
by the difference of the entropy density s(Te) — sh(?e) = (sq(Te) — sh(?e))/2. Although we will focus our arguments 
only on the case of the ansatz ( p!3|) , other forms of interpolation will be necessary when we look into the case with 
different initial conditions for the entropy density. 

Fig. is our final result for the diphoton emission when the system passes through CEP. Because the continuum 
contribution turns out to be 2 '--^ 3 orders of magnitude smaller than the pole contribution, we have plotted only the 
pole contribution with a variety of a ranging from 0.2 (bottom) to 0.9 (top). We cut the invariant-mass M shown in 
the figures up to 170 MeV, for the a mass at the freeze-out temperature is around 170 MeV in CASE-II. 

At a glance we notice a characteristic property that the multiplicity approaches zero as M and has a peak 
around AI = 50 ~ 120 MeV. Although we have expected that the presence of massless a mesons would induce more 
and more significant enhancement in the multiplicity as M gets smaller and smaller, the resulting multiplicity in the 
low energy region is reduced altogether. This suppression is because of the amplitude, or the strength of the pole 
contribution, which is given by Eq. (^. 

The location of the peak actually depends on to* , /* (the a mass and the pion decay constant at (/i = /iE, T = 0)), 
and the evolution of temperature as a function of time. In any case the multiplicity starts from zero at M = 0. As 
the energy gets larger, the amplitude rapidly increases and at the same time the space-time volume in which the 
thermal system stays is reduced by the steeper slope of temperature. Consequently the peak observed in Fig. [t] is 
produced in somewhat general situations due to the behavior of the amplitude and the time evolution. Thus we can 
anticipate that the peak found in our analysis appears in different physical conditions (different choices of to*, /*, 
and the hydrodynamic evolution) as long as the system passes through CEP, though its location is altered according 
to those conditions. 

We would remark that the annihilation process tt+vt^ 27 makes no background in such a low energy region below 
the threshold of two pions. In the scenario of the spectral enhancement in the chiral restoration at zero chemical 
potential, the annihilation process brings about a huge background comparable to the desired signal, which crucially 



smears the spectral changes [18 . In contrast, the main background for the CEP peak comes from the electro-magnetic 
decay into two photons froniTr' in a hot medium and also from tt" emitted after freeze-out. In the present analysis 
the pion mass is fixed at a constant value independently of temperature. In addition to the pole contribution at 
M = TO^o — 135 MeV, the spectral functions in the pion-channel have a continuum contribution in the low energy 
region below the pion mass because of the mass difference between pions and a mesons. However, this pion continuum 
gives rise to so small contribution, which is at most comparable to the continuum contribution in the a-channel, that 
it is sufficiently negligible. 

Although the invariant-mass spectrum has a characteristic shape with a peak whose location is separable from the 
pion pole, the combinatorial background of photons from tt" — > 27 will be crucial for experimental observations. The 
invariant-mass distribution of photon pairs from separate sources is the origin of the combinatorial background. In 
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principle, however, the combinatorial background can be subtracted by means of the event-mixing method. Thus the 
accuracy necessary to identify the peak originating from a mesons is roughly estimated by the comparison between 
the number of photons from a mesons and the number of photons from tt" which might pair with photons from a 
mesons in the construction of the invariant-mass. 

Let us roughly estimate the accuracy necessary to identify the diphoton peak from the a meson in the present 
formulation, though the analysis should not be reliable beyond a qualitative sense and the quantitative results may 
be significantly changed under different conditions for the choice of parameters, the initial values of the temperature 
and the chemical potential, the equation of state, the hydrodynamic evolution, and so on. 

The 7r° yield in the present condition can be estimated by the integration of the Bose distribution function on the 
freeze-out surface as Q 

with Tf = 4.7 fm and Tf — 120 MeV (CASE-I) substituted. This value of the 7r° yield seems small as compared with, 
for example, the particle production observed in AGS |Q. It is because we set the initial point where the system 
reaches equilibrium right at CEP, not above it. As a result the freeze-out time Tf becomes smaller. 

Since almost all the tt^ decays via an electro-magnetic process into two photons, the number of photons from tt*^ is 
estimated as twice of Eq. (|l^). On the other hand, the invariant-mass spectrum of diphoton emission with the range 
up to pt ^ M '--^ 100 MeV integrated out gives 



1 diVcr-y-y 

(Area) AY 



10"^ fm"^ (15) 



As a rough estimate, the accuracy of order ^ 10^ to the 7r° peak after the subtraction of the combinatorial background 
is necessary to detect the clear signal from nearly massless a mesons. The infrared suppression of the amplitude we 
find in the present paper is responsible for the smallness of the resulting output (p^. As we stated above, however, 
the quantitative results here may be considerably changed. For reliable quantitative analyses, we must clarify the 
hydrodynamic properties around CEP, which are just beginning to be investigated p4|. 



V. SUMMARY 



We investigate the spectral functions in the u-channel near CEP where the first-ordered transition of the chiral 
restoration would terminate. The a pole and the continuum from the decay process cr — )■ 2(t dominate over the spectral 
functions near CEP, as expected. Our method has almost no model artifact though it contains some parameters put 
by hand. In the present analysis we try two parameter sets in the extreme cases; one is the case with no medium 
effect and another is the case with a considerably large medium effect. 

Using the resulting spectral functions we have evaluated the multiplicity of diphoton emission. Within the simplest 
space-time evolution described by Bjorken's scaling solution, we present the theoretical prediction for the diphoton 
emission when the system passes through CEP in a hydrodynamic evolution. Our results show that the diphoton 
multiplicity with vanishing transverse momentum has the characteristic shape with a peak around the invariant mass 
50 ^ 120 MeV. We find, contrary to expectation, that the infrared dynamics makes the amplitude of nearly massless 
(J mesons so suppressed that the CEP signal is weakened. The severe background comes from a large peak produced 
by 7r° after freeze-out. 

To proceed further and employ more realistic hydrodynamics, it is essential for the present framework to settle the 
parametrization in the whole (/i, r)-plane and extend it above the critical point. As for the chiral critical point in 
three flavor QCD at finite temperature |^6|, lattice simulations prove to be powerful instruments to search for the 
t-like and the /i-likc directions Unfortunately, however, CEP located at high density is still hard to access by 

means of lattice simulations. 
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APPENDIX A: INTERACTIONS NEAR THE CRITICAL END POINT 

In the presence of finite chemical potential for the baryon density, the effective potential in terms of a and tt can 
be expanded as 

V{a, Tf) = A{a^ + Tr2)3 - ^(a^ + jC^f + C{a^ + _ Da, (Al) 

where the last term embodies explicit breaking of the chiral symmetry due to finite quark masses. The effects arising 
from finite chemical potential induce the first term of the sixth power (finite A) and reduce the second term of the 
fourth power (small B). As a result the phase transition could be first-ordered. In the vicinity of CEP, coefficients A, 
B, C, and D are determined by the following conditions: d^V/dn^ = and d^V/da'^ = for the meson masses 
given by Eq. (|l|), dV/da = for the condensation given by Eq. (^, and d'^V/da^ = for CEP. We need one more 
condition to determine all the coefficients away from CEP because the last condition is peculiar just to CEP. Then 
we assume that A can be regarded as constant in our analysis. Actually we can expect that the value of A is more 
sensitive to /i, which is fixed as /i = /ie throughout this paper, rather than T remembering that A stems from the 
effects of finite chemical potential. Then it follows that 



leading to the full-vertices of the three-point interactions as 



d^V_ ml -ml 



APPENDIX B: FORMULATION OF THE DIPHOTON EMISSION RATE 

The multiplicity of two photons (diphotons) per unit space-time volume is given in the thermal circumstance by 
dTV^ ^ |(/;7(fci,Ai)7(fc2,A2)|5|»)p d^fci d^fc^ 

d'^x ,4"^ Z{l3) ' VT ■ 2u;i(2^)3 ■ 2w2(2^)3' ^ ' 

Al , A2 

where Z{P) is the partition function and the S'-matrix is denoted by S. Two photons have four-momenta (wi = 

(cl'2 = 1^21,^2) and polarization vectors e^(Ai), e^{X2) respectively. The S'-matrix depends on the actual 
processes which can be expressed through the effective Lagrangian, i.e. 

^077 = g<l>fj4>FtJ.iyP'^'^ : (B2) 

with the effective coupling constant g,p-y-y of order ac — 1/137. In our notation cj) denotes collectively any particle 
which can decay into two photons such as 7r° and a mesons. We will focus our attention to the case of ^ = cr in the 
present paper because we are interested in the diphoton spectrum in the low energy region. Up to the lowest order of 
the perturbation with respect to g,p.yy, the ^-matrix can be expanded as 

(/;7(A:i,Ai)7(A;2,A2)|5|z) 

= -4i5^^^{(fci • ■ 6*) - (fci • e^)(fc2 • el)} J d^xe'<^''^+''^Hmix)\t). (B3) 

Since real (on-shell) photons are emitted, we must take the summation with respect to the polarizations only over the 
transverse components. The Ward identity, however, simplifies the manipulations, that is, we reach the correct answer 
simply by taking the summation over all the polarizations. Thus we can make use of the formula, J2\ efiW^KX) — 
—g^iv, to acquire 
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FIG. 8: Diagrams contributing the decay process a ^ 27 in the chiral quark model. The box indicates the full vertex of a ^ 27r 
inferred from the effective potential. 



after taking the summation over states. Here, the Green's function, 
can be related to the spectral function as PSI 



e/3po _ 1 ■ 



(B5) 



(B6) 



where the spectral function is defined as Eq. (|). If the thermal medium has the flow velocity specified by then 
f3pQ in the above expression should be replaced by /3p • u. 

In the kinematics specified by + W2 = w and p = (back-to-back kinematics), we can easily perform the 
integration over ki and ^2. Taking into account the symmetry factor of identical bosons (photons), we finally reach 



dA^ 



(B7) 



APPENDIX C: EFFECTIVE COUPLING CONSTANT 



The effective coupling constant describing the decay a 2j has been calculated within the framework of the linear 
sigma model at zero temperature |]39| and within the framework of the NJL model at finite temperature [^ . In the 
present analysis we estimate the coupling ga-j-y by means of the chiral quark model at finite baryon density as well as 
at finite temperature. Diagrams to be calculated are shown in Fig. ^. The result in the back-to-back kinematics is 
written as 



(CI) 



where ac = e^/47r is the electro- magnetic coupling constant, ti™'^'' and t^'"'"''' represent the vacuum and the medium 
contributions of the pion loops respectively, given by 



, (mat) 



x(mat) 



4mi 



sm 



dq — — '^^-TTT In 



En 



Eq-q 



- 2m^) —nBiuj/2) In 



Ami 



(C2) 
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with the notation Eq = ^ + m^. The vacuum and the medium contributions, tg™'^'' and i^™'**^ come from the quark 
loops, which are given by 



5gmq 



1+ 1 




5.9 771,(^2 — 4m^) 



Eq + q 



Eg-q 



(u) - 2mq) — ^(nF(w/2 + Hq) + nY{uj/2 - fj,q)) In 



LO + ^Ijp- — \m^q 



,2 _ 4^2 



(C3) 



with £^g = y 9^ + TO^ as above, g is the strength of the Yukawa-couphng between mesons and quarks, ruq denotes 

the constituent quark mass. We have chosen g — 3.2 so as to reproduce the constituent quark mass at tree-level, that 
is, mq = gf^ = 300 MeV @. 

It is important to note that diagrams in Fig. ^ have no a loop since the a mesons are electrically neutral. Thus it 
is absolutely necessary to include the contribution from quark loops in contrast to the calculation of the self-energy 
in which the a contribution dominates over the quark loops. 
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